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Abstract 

In this paper we discuss a method for nondestructive measurements of the longi- 
tudinal profile of sub-picosecond electron bunches for X-Ray Free Electron Lasers 
(XFELs) . The method is based on the detection of the Coherent Synchrotron Radi- 
ation (CSR) spectrum produced by a bunch passing a dipole magnet system. This 
work also contains a systematic treatment of synchrotron radiation theory which 
lies at the basis of CSR. Standard theory of synchrotron radiation uses several 
approximations whose applicability limits are often forgotten: here we present a 
systematic discussion about these assumptions. Properties of coherent synchrotron 
radiation from an electron moving along an arc of a circle are then derived and 
discussed. We describe also an effective and practical diagnostic technique based on 
the utilization of an electromagnetic undulator to record the energy of the coher- 
ent radiation pulse into the central cone. This measurement must be repeated many 
times with different undulator resonant frequencies in order to reconstruct the mod- 
ulus of the bunch form-factor. The retrieval of the bunch profile function from these 
data is performed by means of deconvolution techniques: for the present work we 
take advantage of a constrained deconvolution method. We illustrate with numerical 
examples the potential of the proposed method for electron beam diagnostics at the 
TESLA Test Facility (TTF) accelerator. Here we choose, for emphasis, experiments 
aimed at the measure of the strongly non-Gaussian electron bunch profile in the 
TTF femtosecond-mode operation. We demonstrate that a tandem combination of 
a picosecond streak camera and a CSR spectrometer can be used to extract shape 
information from electron bunches with a narrow leading peak and a long tail. 
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1 Introduction 



Electron bunches with very small transverse emittance and high peak cur- 
rent are needed for the operation of XFELs [1,2]. This is achieved using a 
two-step strategy: first generate beams with small transverse emittance using 
an RF photocathode and, second, apply longitudinal compression at high en- 
ergy using a magnetic chicane. The bunch length for XFEL applications is of 
order of 100 femtoseconds. Since detailed understanding of longitudinal dy- 
namics in this new domain of accelerator physics is of paramount importance 
for FEL performance, experiments on this subject are planned in test facili- 
ties. The femtosecond time scale is beyond the range of standard electronic 
display instrumentation and the development of nondestructive methods for 
the measurement of the longitudinal beam current distribution in such short 
bunches is undoubtedly a challenging problem: in this paper we discuss one 
of these methods, which is based on spectral measurements of Coherent Syn- 
chrotron Radiation (CSR) produced by a bunch passing a dipole magnet or 
an undulator. 

After the compression stage, the CSR pulse can be analyzed, for example, 
with a spectrometer. As we point out in Section 2, we can decompose the 
coherent radiation spectrum, P(u), into the product of the square modulus of 
the bunch form factor, | F(u>) | 2 , and the single particle radiation spectrum, 
p(ou) [3]. The CSR spectrum, then, provides information on the bunch form 
factor, although one has to keep in mind that, in order to find | F(u) | 2 from 
P{oj), one needs the quantity p(cu). 

Section 3 contains a systematic treatment of synchrotron radiation theory. 
All the results presented in this section are derived from fundamental laws of 
electrodynamics and the reader can follow the whole derivation process from 
beginning to end. We use a synthetic approach to present the material: simple 
situations are studied first, and more complicated ones are introduced gradu- 
ally. Previous experience in synchrotron radiation theory would be helpful but 
is not absolutely necessary, because all the required material is independently 
derived. In this respect our paper is reasonably self-contained. 

A clear message in this work is that reexamination of dogmatic "truths" 
can sometimes yield surprises. For years we were led to believe that famous 
Schwinger's formulas [4] are directly applicable to the case of synchrotron ra- 
diation from dipole magnet and even now no attention is usually paid to the 
region of applicability of these expressions. While such formulas are valid in 
order to describe radiation from a dipole in the X-ray range, their long wave- 
length asymptotic are not valid, in general. Analytical study of this matter 
was first performed in the 80's [5]. However, standard texts on synchrotron 
radiation theory do not seem to provide a derivation of the expression for the 
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synchrotron radiation from an electron moving along an arc of a circle. In 
Section 3 we present such a derivation, which we believe is quite simple and 
instructive. Properties of coherent synchrotron radiation from dipole magnets 
in the time domain are derived and discussed too. 

Standard theory of synchrotron radiation relies upon other approximations 
too, and it seems interesting to pay attention to their region of applicability. 
To be specific, two important limitations are discussed. First, it is usually 
assumed that the observer lies at infinite distance from the source. Second, 
people don't pay attention to the fact that, in real experimental conditions, 
the radiation is seen by the detector through some limited aperture. Thus, 
both finite distance effects and diffraction effects are ignored. The aperture 
sizes and the distances for which the latter assumptions are valid are so large, 
in long wavelength range, as to be of limited practical interest. We can draw 
at least two main conclusions from Section 3: first, long wavelength radiation 
spectrum distortions can arise, physically, from violation of the far zone as- 
sumption or from aperture limitations (or both). Second, a CSR diagnostic 
method based on the famous Schwinger's formulas (or the ones derived in 
[5]) is of purely theoretical interest. In real accelerators, the long wavelength 
synchrotron radiation from bending magnet in the near zone integrates over 
many different vacuum chamber pieces with widely varying aperture and p(u) 
is usually a very difficult quantity to calculate with great accuracy. 

An effective and practical technique based on the spectral properties of 
undulator radiation can be used to characterize the bunch profile function. 
The method we describe in Section 4 uses an electromagnetic undulator and 
it is based on recording the energy of coherent radiation pulses in the central 
cone. This coherent radiation energy turns out to be proportional, per pulse, 
to the square modulus of the bunch form-factor at the resonant frequency of 
the fundamental harmonic. The measurement must be repeated many times 
with different undulator resonant frequencies (which are tuned by changing 
the undulator parameters) in order to reconstruct the modulus of the bunch 
form-factor. 

The retrieval of the bunch profile function from the modulus of its form 
factor is preformed, in Section 5, by means of a deconvolution technique. 
For the present work we choose a constrained deconvolution method. This 
consists in finding the best estimate of the bunch profile function, F(t), for 
a particularly measured form-factor modulus, | F{uj) |, including utilization 
of any a priori available information about F(t). In the end of the paper we 
illustrate with a numerical example the potential of the proposed technique 
for diagnostics at the TESLA Test Facility accelerator. Here we have chosen, 
for emphasis, an experiment aimed at measuring the strongly non-Gaussian 
electron bunch profile at TTF, Phase 2 in femtosecond mode operation. The 
femtosecond mode operation is based on the experience obtained during the 
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operation of the TTF FEL, Phase 1 [6] and it requires one bunch compressor 
only. An electron bunch with a sharp spike at the head is prepared with an 
rms width of about 20 microns and a peak current of about one kA. This 
spike in the bunch generates FEL pulses with duration below one hundred 
femtoseconds. We demonstrate that a tandem combination of a picosecond 
streak camera and CSR spectrometer can be used to extract shape information 
from electron bunches with a narrow leading peak and a long tail. 



2 Physics of coherent synchrotron radiation 

To begin our consideration, let us recall some well-known aspects of CSR. 
From a microscopic viewpoint, the electron beam current at the entrance of 
a bending magnet is made up of moving electrons arriving randomly at the 
entrance of the bending magnet: 



where S(-) is the delta function, (— e) is the (negative) electron charge, N 
is the number of electrons in a bunch and tk is the random arrival time of 
the electron at the bending magnet entrance. The electron bunch profile is 
described by the profile function F(t). The beam current averaged over an 
ensemble of bunches can be written in the form: 



The profile function for an electron beam with Gaussian current distribution 
is given by: 



and the probability of arrival of an electron during the time interval (t,t + dt) 
is just equal to F(t) dt. 

The electron beam current, J if), and its Fourier transform, J(u>), are con- 
nected by 



N 



•/(*) = (-e) , 



fc=i 



(J(t)) = (-e)NF(t) . 




oo 




— oo 



exp(iut)J(t)dt = (-e) ^2 exp(icjt fc ) , 

fc=i 
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J(t) = — [ &v(-icjt)J(u)dw = (-e)Y,6(t-t k ) 
271 J k=i 



and, therefore, the average value of | J{oS) | 2 can be written as: 
{J{uo)J*{uo)) = e 2 AT + e 2 ^(exp(i^ fc ))(exp(-icc;t„)) . 



The expression (exp(io;^)) is nothing but the Fourier transform of the bunch 
profile function F(t), since: 

oo 

(exp(iu;t k )) = J F(t k ) exp([Lut k ) dt k = F(lu) . 

— oo 

Thus we can write: 

(| J(uj) | 2 ) = e 2 N + e 2 N(N - 1) | F(u) | 2 , 

where the Fourier transform of the Gaussian profile function has the form: 
FH = exp (-^f) . (1) 



Above we described the properties of the input signal in the frequency 
domain. The next step is the derivation of the spectral function connecting 
the Fourier amplitudes of the output field and the Fourier amplitudes of the 
input signal. We will investigate the synchrotron radiation in the framework 
of a one-dimensional model. Also, we will assume that the cross section of the 
particle beam is small compared with the distance to the observer, so that 
the path length difference from any point of the beam cross section to the 
observer are small compared to the shortest wavelength involved. 

The component of synchrotron radiation electric field in time domain, E Xjy (t), 
and its Fourier transform, E Xty (u), are connected by 



1 °° 

= 7^ J E XtV (u;)ex.p(-iu;t)du> , 



and the Fourier harmonic for oo < is defined by the relation E*(uj) = E(—uj). 
On the other hand, the Fourier harmonic of the electromagnetic field and the 
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Fourier harmonic of the current at the dipole magnet entrance are connected 
by: 

E XtV (u) = A Xty (uj)J{uj) , lu > , 

where A Xty (u) is the synchrotron radiation spectral function of the dipole. 
Since the radiation power is proportional to the square of the radiation field, 
the averaged total power at a certain frequency u> is given by the expression: 

P(u) = p(u) [N + N(N -1)| F(u) | 2 ] , (2) 

where p(u) is the radiation power from one electron. The first term N in 
square bracket represents ordinary incoherent synchrotron radiation with a 
power proportional to the number of radiating particles. The second term rep- 
resents coherent synchrotron radiation. The actual coherent radiation power 
spectrum depends on the particular particle distribution in the bunch. For 
photon wavelengths equal and longer than the bunch length, we expect all 
particles within a bunch to radiate coherently and the intensity to be pro- 
portional to the square of the number of iV of particles rather than linearly 
proportional to N, as in the usual incoherent case. This quadratic effect can 
greatly enhance the radiation since the bunch population can be from 10 8 to 
10 11 electrons. On the other hand the coherent radiation power falls off rapidly 
for wavelengths as short or even shorter than the rms bunch length. 

A method for estimating the bunch profile function F(t) is to compute the 
factor p(u), to measure the CSR power spectrum P(u) and to deduce the 
factor | F(u) | 2 from equation (2). However, it is clear that the retrievable 
information about the profile function will, in general, not be complete, for 
it is the particular measured square modulus of the form-factor that can be 
obtained, not the complex form-factor itself. 



3 Radiation from dipole magnet 

The phenomenon of coherent synchrotron radiation has been introduced in 
a conceptual way in the preceding section. As was pointed out, the CSR power 
spectrum is the product of the square modulus of bunch form-factor | F(cu) | 2 , 
and the single-particle power spectrum p(cu). Any diagnostics method based on 
CSR must devote therefore attention to correctly determinate the synchrotron 
radiation spectrum from one electron p(ou), for this quantity plays a critical 
role in the form-factor measurement. In the present section we will deal with 
the properties of the single-particle radiation in time and in frequency domain 
as well as with the CSR properties in time domain. In order to give the sub- 



5 



Fig. 1. Geometry for synchrotron radiation production from a bending magnet 



ject a semblance of continuity, it will be desirable to introduce considerable 
matter which can be found in any of standard text on synchrotron radiation 
theory (see for example [3] , [7] , [8] ) . The typical textbook treatment consists in 
finding the expression for the synchrotron radiation spectrum from an electron 
moving in a circle. However, no attention is usually paid to the region of ap- 
plicability of the derived expressions. For instance, the standard extension of 
the theory to the case of synchrotron radiation from dipole magnet is based on 
the assumption that the energy spectrum formula is equivalent to the famous 
Schwinger's formula [4]. While this formula is valid for the X-ray range, it 
does not provide a satisfactory description in the long wavelength asymptotic 
[5]. Moreover Scwinger's formulas are found in the limit of an infinite observer 
distance and under the approximation of no limiting aperture through which 
the radiation is collected. Here we pay particular attention to the region of 
applicability of Schwinger's formulas; specifically, we derive expressions for the 
time and frequency dependence of the electromagnetic radiation produced by 
an electron moving along an arc of a circle, and we investigate how the CSR 
field pulse in the time domain is modified in more realistic situations. 



3. 1 Radiation field in the time domain 

Let us start reviewing how a single electron radiates as it moves in a dipole 

magnet. The electron is subjected to an acceleration, v = c/3, of magnitude 
c 2 f3 2 / R so that a distant observer detects the electric and magnetic fields in 
the form of electromagnetic waves from the radiating electron. Fig. 1 shows 
the relationship between the observer at a fixed point P, whose coordinates are 
(r*o, t), and the radiating electron at (R, t'), t' being the emission (or retarded) 
time. The fundamental laws of electrodynamics tell that the electric field of 
a charge (— e) moving along an arbitrary trajectory is given by the Lienard- 
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Fig. 2. Synchrotron radiation production from circular motion. Directions of accel- 
eration vector as seen by an observer at P, drawn in plane view 




Fig. 3. Time variation of a synchrotron radiation pulse generated by a highly rela- 
tivistic electron moving in a circle as seen by an observer in the orbital plane 

Wiechert formula 



E(t) = 



-e)J_ 



+ 



(-e) 



c r 



n x [(ft - (3) x (3} 
(1 - n-/3) 3 



where n is a unit vector along the line from the point at which the radiation is 
emitted at the emission time to the observation point at the observation time, 
and we understand that the quantity in brackets must be evaluated at the 
retarded time i! = t — - c \ f(t') | . The latter equation consists of two distinct 
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R0 



Fig. 4. Synchrotron radiation production from an arc of a circle. Directions of ac- 
celeration vector as seen by an observer at P, drawn in plane view 



E x (A.U.) 




-0.2-1 



Fig. 5. Time variation of a synchrotron radiation pulse generated by a highly rel- 
ativistic electron moving along an arc of a circle. The normalized bending angle is 

4>m = 70m = 2.5 

parts. The first is inversely proportional to the square of the distance between 
radiation source and observer, depends only on the charge velocity and is 
known as velocity or Coulomb field. The second is inversely proportional to 
the distance from the charge, depends also on the charge acceleration and it 
is known as acceleration or radiation field. At large distances from the moving 
electron, the acceleration-related term dominates, and it is usually associated 
to the electromagnetic radiation of the charge. The region of space where 
the radiation field dominates is called a far (or wave) zone and the radiative 
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E (A.U.) 




Fig. 6. Time variation of a synchrotron radiation pulse generated by a highly rel- 
ativistic electron moving along an arc of a circle. The normalized bending angle is 

<t>m = 70m = 1-5 

electric field in the far zone is given by the formula 



c r 



nx[(n-(3)x 0] 
(1 -n-(3f 



(3) 



We can use (3) to look at all kinds of interesting problems. This is a compli- 
cated expression, but it is easy enough to be used in a computer calculation 
which can be further visualized as a geometrical picture. Such geometrical 
pictures will give us a good qualitative description of the situation. Usual the- 
ory of synchrotron radiation is based on the assumption that the electron is 
moving on a circle and radiation is observed from the whole circular trajec- 
tory, so that in each cycle we get a sharp pulse of electric field. A far-field 
computation of the predicted time dependence of synchrotron radiation for 
circular motion is presented in Fig. 3. The horizontal component of electric 
field plotted versus the normalized variable uj c t, where uo c = 37 3 c/ (2R) is the 
critical frequency of synchrotron radiation. The field in the orbital plane has 
a zero around t = a;" 1 . Numerically from Fig. 3 one obtains 

oo 

J E r dt = . (4) 

— oo 
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Magnet Length = 0.5 



o.o 



-0.2 



cot 



Fig. 7. Time variation of a synchrotron radiation pulse from a short bending magnet. 
The normalized bending angle is 4> m = 7^ m = 0.5 



observer 
P 




uniform current 



Fig. 8. Geometry for an infinite uniform current progressing through an arc of a 
circle and relative position of an observer 

It is interesting to stress the fact that equation (4) is strictly related to the well- 
known result that a uniformly charged ring does not radiate. In fact, starting 
from (4) one can show that a system of N identical equidistant charges (— e) 
moving with constant velocity v along a circle does not radiate in the limit 
for N — > oo and (—e)N = const, and the electric and the magnetic fields of 
the system are the usual static values. 

It is important to realize that (4) is valid only when the electron is moving 
in a circle. But we now want to study synchrotron radiation from an electron 
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moving along the arc of a circle (see Fig. 4). By inspecting Fig. 5-7, one 
can see that when the electron moves in arcs of circles with different angular 
extensions <p m , the time-average of the electric field is nonzero. 

Let us consider for a moment (in parallel with what has been done in the 
case of a circle) the case of infinitely long electron bunch with the homoge- 
neous linear density Ao and current / = — ev Xq- A current circuit consists of 
the arc and semi-infinite straight lines (see Fig. 8). The angle between the 
straight lines is equal to the bending angle m of the magnet. The fact that 
the average electric field from a single particle is different from zero means 
that the acceleration field from our infinite circuit must be different from zero 
too. Then we have deal with an intriguing paradox, since it is a well known 
result that an uniform electron current does not radiate, not only in the case 
of circular motion, but independently from the trajectory 1 . 

It is possible to explain this contradiction in very simple terms as follows. 
First of all we know that the velocity (electric) field from a line current (in- 
cluding our case) is proportional to 1/ ] f± |, where | r_j_ | is the distance of 
the observer from the line charge. Second, in our case, the acceleration part 
of the electric field is proportional to 1/ | f |, where | r | is the distance of the 
observer from the magnet (since the acceleration field sources are strictly lim- 
ited to the particles in the magnet only). The situation is depicted in Fig. 8: 
the ratio | f± \ / \ r \ is finite for any position of observer therefore there is no 
region in space where the acceleration field dominates. This last observation 
solves our paradox since, in the case of infinite charge current in an arc of a 
circle, we cannot talk about far zone at all, although a non-zero acceleration 
field is present. 

It is interesting to discuss, at this point, the shape of the coherent syn- 
chrotron radiation field pulse under different trajectories. What we have been 
dealing with before is simply a far field analysis of the single-particle radiation 
field. When a large number of electrons move together, all the same way, the 
total field will be a linear superposition of the individual particle fields. Of 
course the result depends upon the longitudinal distribution of the electrons. 
One aspect of the problem that we can immediately deal with is the coher- 
ent synchrotron radiation production from a "short" magnet, i.e. <p m < I/7. 
We can assert that the field appears as shown in Fig. 9. In fact the time- 
dependence of the field has, for every electron, the shape shown in Fig. 7 and 
the total field emitted by the electron bunch is represented by a sum of these 
pulses, one for each radiating electron. In the limit for "short" magnets, the 
reader can easily come to the intuitive conclusion that the time profile of the 



It can be proved that the radiative interaction force in the longitudinal direction 
(parallel, at any time, to the velocity vector by definition) is equal to zero at any 
point of such circuit 
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Fig. 9. Coherent synchrotron radiation production from a "short" magnet. The 
time profile of the electron bunch density is linearly encoded onto the electric field 
of the radiation pulse. Spikes correspond to spontaneous radiation, and average 
value of electric field corresponds to CSR. The width of the temporal profile of the 
electric field corresponds directly to the electron bunch length, and the shape of the 
temporal profile is proportional to the longitudinal distribution within the bunch 

electron bunch density is linearly encoded onto the electric field of the radia- 
tion pulse. The width of the temporal profile of the electric field corresponds 
directly to the electron bunch length, and the shape of the temporal profile is 
proportional to the longitudinal bunch distribution. 

We have just argued that the bunch density is linear encoded onto the elec- 
tric field as in Fig. 9; nevertheless, in Fig. 9, we illustrate small fluctuations 
which occur in the field amplitude too. The reason for this is that the electron 
bunch is composed of large number of electrons, thus fluctuations always exist 
in the electron beam density due to shot noise effects. For any synchrotron 
radiation beam there is some characteristic time, which determines the time 
scale of the random field fluctuations. This characteristic time is called co- 
herence time r c of synchrotron radiation and its magnitude is of the order 
of the pulse duration from one electron. The physical significance of these 
fluctuations is that there is a short wavelength radiation component of the 
radiation in the range of the inverse pulse duration from a single electron. 
Simple physical considerations show that the energy spectrum of this hard 
radiation component is order of P(uj) — Np(u). This explains the relation 
between small fluctuations of the radiation field amplitude and spontaneous 
(incoherent) emission of synchrotron radiation. 

Up to this point we only talked about coherent radiation from a "short" 
magnet. Ultimately we want to consider CSR production from an arc of a 
circle. An easier step in this direction consists in the analysis of the CSR time 
pulse from a circular motion. It is possible, indeed, to derive a simple analytical 
expression for the CSR pulse from a bunch with an arbitrary distribution of 
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the linear density satisfying the following condition: 
R dF(t) 



C7 3 dt 



« F(t) . (5) 



The latter condition simply indicates that the characteristic length of the 
bunch is much larger than -R/7 3 . Let us express the total CSR pulse as a 
superposition of single particle fields at a given position in the far zone: 

00 

£csr(*)= / E r {t-r)NF{r)AT , (6) 

—00 



where we calibrated the observer time in such a way that, when F(r) = S(t), 
the single particle radiation pulse has its maximum at t = 0. To calculate 
the integral in (6) one should take into account the property (4) of the kernel 

— * 

E T (t — t) which has been discussed above. Using (4) and (5) one can simplify 
equation (6) in the following way. The integral (6) is written down as a sum 
of three integrals 

t-Si 

Ecsn{t)= J E r (t-r)NF(r)dr 

—00 

t+&2 OO 

+ J E r (t — t)NF(t) d t + J E r (t-r)NF(r)dr , (7) 

t-Si t+8 2 



where 5( ij2 ) satisfy the following conditions: 



*(i,2)»^, S {1 , 2) ^ « F(t) . (8) 



The bunch profile function is, therefore, a slowly varying function of the time 
and we simply take F(t) outside the integral sign and call it F(t) when cal- 
culating the second integral of (7): 

t+S-2 t+6 2 

J E T (t — t)NF(t) (It - NF(t) J E r (t - r) dr . (9) 

t-Si t-Si 



Then we remember that the average of the electric field over time is zero when 
an electron is moving in a circle. As a result, we rewrite the integral in the 
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form: 



t+s 2 



J m-r) 



dr = - 



t-Si 



t-Si 



J E T {t-r)dT+ J E T {t-r)dr 

-oo t+82 



(10) 



Taking into account (9) and (10), the expression for CSR pulse becomes 



t-Si t-Si 

E C sn(t)= I E r (t-r)NF(r)dr-NF(t) f E T {t - r) dr 



+ / E r (t-T)NF(r)dT- NF(t) J E T (t - r) dr 

t+82 t+82 



'ID 



Integrating by parts, the first pair of integrals on the right hand side of (11) 
can be joined in a single one; the same can be done with the second pair: 



t-81 



r r-i dF(r) 
£csr(*) = - / * E r ( t -r)N^—^d 
J dr 



Er 



t + 82 



dr 



;i2) 



where $ E r is a primitive of E T . What is left to do now is to evaluate a 
primitive of the radiation field from one electron. 

To calculate the primitive of E T we use (3) and note that the electric field 
is expressed in terms of quantities at the retarded time t'. The calculation 
is simplified if we use the following consideration: since, in general dt/ dt' = 
(1 — n ■ j3) and 



n x [(ft — 13) x f3] 
(1 - n-/3) 2 



d 


n x [(nx (3)\ 




dt' 


(1-ft-P) 




we also have 




d 


'fix [{fix (3)] 




dt 


(1-n-P) 





nx[(n-/3)x f3] 

' (i-fi-pf 



Thus we can write -Bcsr(^) as 
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Fig. 10. Geometry for synchrotron radiation from circular motion 



^csr(0 = 



+ 



c r 



t+&2 



t-Si 



r 



n x [(nx (3)} 
(l-n-P) 



J (t-T) 



dr 



n x [(n x /?)] 



{t-T) 



dr 



(13) 



where the quantity in brackets must be evaluated at the retarded time i! = 

(t-T)-±\f(t')\. 

Now we would like to compute the quantities required for (13). Since (8) 
holds, we can substitute the function in brackets in both integrals on the right 
hand side of (13) with its asymptotic behavior at r > R/(cy 3 ). Because the 
angles are very small and the relativistic 7 factor is very large, it is very useful 
to express (13) in a small angle approximation. The triple vector product is 
calculated from Fig. 10 



n x [n x (3\ = u t'e x + 



(14) 



Here 6 is the vertical angle, uj = c/R is the revolution frequency and e Xty are 
unit vectors directed along the x and y axis of the fixed Cartesian coordinate 
system (x, y, z) shown in Fig. 10. The definition of n and R can be used to 
compute the scalar product in the denominator in (13) so that 



n ■ p = p cos 6 cos LU t' ~ (3(1 - 9 2 /2) (1 - (cu t') 2 /2) . 



(15) 



We assume, here, that the vertical angle is small enough and we can leave out 
the cosine factor. We can now write E CSR (t) as 
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£csr(0 



2e 



t-5i r 



+ 



t+8 2 



c | r 

(^0^4 + 9e y 
MO 2 



(u t'e x + 9e y 



(t-T) 



dr 



(t-T) 



dr 



(16) 



Yet, part of the integrand in (16) is still expressed as a function of t', which 
has to be converted in a function of t — r using the explicit dependence 

£ — r = £' H — I r(t') 1= t' cos 9 sin cu t' . 

c c 



Since we assume the vertical angle is very small, we may use the replacement 
cos 9 ~ 1 . We can therefore approximate t — t by 



1 



t' + - | r(0 | = 

c c 



fo| +ni-/3) + "° W 



6 



We conventionally fixed E v (0) as the maximum value of the field (in time) and 
we are interested in the asymptotic for r ^> Rj (cy 3 ) only, therefore we can 
simply write t — r = cjQ(t') 3 /6 . Solution of this equation allows us to write 
equation (16) as 



-^csr(^) = 



2eN 
c I r I 



[6u (t - r)] 1 / 3 ^ + 9e y \ dj p( r 



(t+82 



[6a; (r-t)] 2 / 3 



dr 



dr 



t-Si 



+ 



[Quoit -T)] 1 ^e x + 9e y ] dF (r 



[Qu {t -r)] 2 /3 



dr 



dr 



(17) 



Limitation (8) indicates that the contribution from the integrands in the right 
hand side of (19) are negligible in the region (t + 62, t — 5\) and therefore we 
can rewrite (19) in its final form: 



EcsR(t) 



2eN 
c I r 



r e{t — r)e x 9e y 

J [6U I t-T |]V3 + [Q Uq I t-T |] 2 / 3 
-00 



dF(r) 
dr 



dr .(18) 



where 

e(t - t) = 1 for < (t - t) < 00 , 
e(t-T) = -l for - 00 < (t - t) < 
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, E x( A -U-) 




Fig. 11. Time structure of a CSR pulse from a Gaussian electron bunch moving in 
a circle. Here 8 = 0. Circles present the results obtained from direct superposition 
of single particle pulses, while the solid line corresponds to the shape calculated by 
means of (18). The rms bunch length is <tt = 100A c /c, where A c = 47ri?/(3c7 3 ) 



It might be worth to remark that the ultimate reason for using the auxiliary 
times 5(i,2) m the derivation of (18) is that they help recognizing the validity 
of the asymptotic substitution, since otherwise a direct integration by parts 
of (6) would immediately give 



oo 

-^csr(^) = I -, I / 
c\r \ J 

— oo 



Under the accepted limitation on the axial gradient of the beam current, this 
equation transforms to (18). 

As an example we show how to use (18) in order to calculate the CSR 
pulse. Let us concentrate on the CSR radiation produced in the orbital plane. 
In this case = and it is obvious that the radiation for such an observer 
is horizontally polarized. To be specific, we consider an electron beam with a 
Gaussian axial profile of the current density 

F(()= 7ik exp hM • 
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where ctt is the rms electron pulse duration. The rms is assumed to be large, 
o" T ^> i?/(c7 3 ). When 9 = and bunch profile is Gaussian profile, we can 
write (18) in the form 



2(-e)N 



oo 




Hi) dT - ^ 



E x (t) 



(27T) 1 /2 6 l/3 CT 3 tu l/3 c | ^ 



/ 



— oo 



Fig. 11 presents comparative results obtained by means of analytical calcula- 
tions (solid curve, calculated with (19)) and numerical results (circles, calcu- 
lated by direct superposition of single pulses from (6)). It is seen from these 
plots that there is a very good agreement between numerical and analytical 
results. 

It is relevant to make some remarks about the region of applicability of 
(18). It is important to realize that (18) is valid only when the electrons 
are moving in a circle and the observer is located in such a way that both 
the velocity term in Lienard-Wiechert formula can be neglected and the unit 
vector n can be considered constant. Another basic assumption is that the 
current density changes slowly on the scale of -R/7 3 . As a rule, this condition 
is well satisfied in all practical problems. It should be also mentioned that 
the above expressions are good approximations only for small enough vertical 
angles (even though they may be immediately generalized). In fact we used 
the replacement cos 6 ~ 1 in the retardation equation and, in practice, such an 
assumption is valid for the range 9 2 <C (cr/i?) 2 / 3 , where a is the characteristic 
length of electron bunch. 

Now we want to extend the latter results to the case of an arc of a circle. 
At this point we find it convenient to impose the following restriction 2 : we 
focus only on the radiation seen by an observer lying at large distance from 
the sources, on the tangent to the electrons orbit at the middle point of the 
magnet. In this case we can continue to use the fixed coordinate system (x, y, z) 
shown in Fig. 10. The observation point and the vector n are within the (y, z)- 
plane and radiation is emitted at an angle 6 with respect to the z-axis. Let us 
start expressing the total CSR pulse as a superposition of single particle fields 
at the given position in the far zone. In the case of an arc of a circle equation 
(6) modifies as follows: 



2 The reader can wonder why it is necessary to describe this particular situation. 
The answer is that this is a difficult subject, and the best way to study it is to do it 
slowly. Although here we deal with a particular example, all the expressions which 
we derive are immediately generalizable 



t+T 




(20) 
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Here the time T in the integration limits is in loco of a window function in 
the integrand, in order to cut the contributions of the single particle radia- 
tion pulse when the electron is not in the arc. This expression contains the 
observation time T, which should be replaced by the retarded time t' e . The 
two times are related by 

2T^t' e + -\f(t' e )\--\f \, 
c c 

where t' e = m /(c<Jo) and m is the bending magnet angular extension. Our 
analysis focuses on the case of a long bending magnet, 70 m ^> 1. Using (4) 
and (5), the field of the CSR pulse is readily shown to be 



t-s 



E C sn(t)= J E T (t — t)NF(t) d t — NF(t) J E r (t-r)dr 

t-T -oo 
t+T oo 

+ J E T (t - t)NF(t) d t - NF(t) J E T {t-r)dr. (21) 



As we have already done previously, we assume that 5 satisfies condition (8). 
Adding and subtracting suitable edge terms one can still perform integration 
by parts, thus obtaining: 



t+T t-T 

E CSK (t)=NF(t + T) J E T (t-r)dr-NF(t-T) J E r (t - r) dr 

— CO —CO 

t-S „, s t+T 



- J (t-T)N^-jp-dT- J $ 

t-T '' " t+S 



dF_(r) 
d' 



E r 



dF(r) 

(t-T)N^-^dr . (22) 
dr 



Since condition (8) holds for 5 we may substitute the 3rd and the 4th integral 
in (22) with a single integral in which the primitive, $ E T , is substituted by its 

asymptotic for large values of the argument, $ E^ . Under the assumption 
of a long bending magnet (uj c T ^> 1) the 1st and the 2nd integral can be 
expressed by means of the primitive asymptotic too. Moreover, we can perform 
a change of variables in all the integrals t — r — > r. As a result expression (22) 
can be written in the form: 



T 



E CS R(t) = —NF(t + T) J 4 A (r) d r - NF(t - T) J E^(r) d ' 

—T oo 
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-T 



(23) 



Using the ultrarelativistic approximation we can calculate a primitive $ 
using (14) and (15). Again we assume that the vertical angle is very sma' 
that we may use the replacement cos# ~ 1. In this situation we have 



E r 
1 so 



n x [(nx P)] iu t'e x + 9e y 

This quantity must be evaluated at the retarded time t' ~ [Qt/uJq] 1 ^. Substi- 
tution these expressions in (23) gives 



^csr(0 



2eN 
c | r 



£{T)e, r 



+ 



[6u | r |]V3 [q Uq | r |] 2 / 3 



dF(t-r) 
dr 



dr 



+ [F(t + T) + F(t - T)] 
- [F(t + T) — F(t - T)] 



(6c T)V3 
(6u T)^ 



(24) 



where T = <^/(12w ). 

As an example of the application of this expression, consider the situation 
when 9 = and bunch profile is a Gaussian. According to (24) the CSR field 
in this case is given by 

E x (t)=G 1 (t) + G 2 (t) 



where 
G 1 

Go 



2(-e)N 



( 27r ) l/2 6 l/3o-3 ^/3 C | f() 

2eiV 

(27r) 1 /2 6 V3 (7T ( a ; T)i/3 c I f Q 



e(r)(t-r) / (t - r) 2 \ 

; exp - v „ ; dr, 



1 1/3 



2a 2 . 



/ (t + T) 2 ^ / (t-T) 2 ^ 

exp - 2 + exp 



2a 2 



2a£ 



Fig. 12 presents the results of calculations with the analytical formula (24). 
Comparison with numerical results shows a very good agreement. 
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Fig. 12. Time structure of the CSR pulse from a Gaussian electron bunch moving 
along an arc of a circle. Here 6 = 0, p = <^/(6u;o0"t) = 2. The physical meaning 
of the magnet length parameter p will be clear after reading Section 3.3. Circles 
present the results obtained from direct superposition of single particle pulses. The 
solid line corresponds to the shape calculated by means of (24) 

3.2 Radiation field in the frequency domain 



Let us now go back to our single particle treatment and proceed to the 
calculation of the spectrum of the electromagnetic radiation produced by an 
electron during a dipole magnet pass. We know that, in general, for any kind 
of motion 

d 2 W _ c\f\ 2 \ E(t) | 2 

d ft d t ~ 4:71 



Here | E | 2 is the squared modulus of the electromagnetic field vector at 
the observation point. The stationary observer who is detecting the radiation 
emitted into the solid angle dft measures the total energy as 



d W c\r 



dft 4vr 



/ I E T (t) \ 2 dt , 



where we used the fact that, in the far field approximation, the distance | f \ 
is much larger than | R | , and the zero order expansion reads | f(t') | = | r | 
(see Fig. 1). From the discussion above, we know that the radiation pulse is 
emitted over a very short period of time so that the only finite contribution 
to this integral comes from times close to t — 0. Extension of the integral to 
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infinite times is only a mathematical convenience which does not affect the 
physical result. 



The general method to derive the frequency spectrum is to transform the 
electric field from the time domain to the frequency domain. Expressing the 
electrical field E T (t) by its Fourier transform, we set 



E r (uj) = J E r (t)exp(-iujt)dt 



(25) 



Applying Parseval's theorem we have 

oo 

E T {u) | 2 dw = 2tt f | E T (t) | 2 d* 



and the total absorbed radiation energy from a single pass is therefore 
dW 



dtt 8tt 2 



C I - 12 

r 



J \Et(u) 



duj . 



Evaluating the electric field by Fourier components we derive an expression 
for the spectral distribution of the radiation energy 



d 2 W c 



ro | 2 | EM | 2 



(26) 



In order to calculate the Fourier transform we can use (3) and note that 
the electric field is expressed in terms of quantities at the retarded time t' . 
The calculation is simplified if we express the whole integrand in (25) at the 
retarded t' remembering t' — t — - \ r(t') |: 



EM = 



c r 



n x [(n — 13) x j3] 
(l-n-/3) 2 



exp[icj(t'+ | f | /c)\dt' 



— * 

Since the distance | f{t') \ is much larger than | R \ , the zero order approxi- 
mation would make | r{t') \ = \ r | (see Fig. 1). However, such approximation 
is not good enough in the exponential factor and we must take the next ap- 
proximation order writing | r(t') \ = \ tq | — n ■ R(t'). Therefore we get 

oo 

c r J 

— oo 



n x [(ft — (3) x (3} 
(l-n-/3) 2 



exp[iu(t f - n ■ R{t')/c)\ dt' . (27) 
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We established some basic equations in the frequency domain with which we 
can start calculations for individual problems. Analytical calculations can be 
performed without big difficulty in two limiting cases, namely the case of 
circular motion and "short" magnet. 



3.2.1 Radiation from circular motion 



In his paper on synchrotron radiation, Schwinger gives remarkable formulas 
for the radiation spectrum in the case of an electron moving in a circle [4]. 
One can find textbooks telling that Schwinger's formulas apply to the analysis 
of synchrotron radiation from an electron moving along the arc of a circle too 
(see for example [8]). This extension is not a physical law: it is merely the 
statement of an approximation which is valid about the entire wavelength 
range interesting for the ordinary user. In this context, a critical study of 
the theoretical status of Schwinger's formulas seems to be of considerable 
importance. 

We are going to apply (27) to our analysis of synchrotron radiation from 
circular motion taking advantage of an expression we previously used in Sec- 
tion 3.1: in the case of circular motion E T (u) can be evaluated remembering 
that 



_d_ 
dt' 



n x [(nx (3)] 



nx[(fj-(3)x 0] 
(1 - n-/3) 2 



Integration by parts yields 



fx x [(n — (3) x (3] 



l-n-(3) 2 

— * 

fix [fix f3] 

■ exp 



exp[iw(t' -nR(t')/c)]dt' 

t'=oo 

exD \\uj(t' -n- R(t'Mc\\ 

(l-fi-(3) 

oo 

ito / n x [n x /3] exp \uj(t' — n ■ R(t')/c) dt' . 



t'=-oo 



(28) 



The reason for using integration by parts is that the contribution from the 
first term (which to be evaluated at t! = ±oo) is zero. Thus we can write 



d 2 W e 2 uj 2 

d^d^ 47T 2 C 



J nx [fix (3] exp [iu;(t' - n ■ R(t')/c)\ dt' 



(29) 
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We must emphasize, however, that this expression is not valid in general, but 
only when a particle is moving in a circle. 



The integrand in (29) can be expressed in components to simplify the inte- 
gration. If we look in the plane of the circle, radiation does not depend upon 
the azimuthal angle. There is a better way to write out the integrand in (29) 
by making use of the azimuthal symmetry, which allow one to use again the 
fixed coordinate system (x,y,z) shown in Fig. 10. The observation point is 
far away from the source point and we focus on the radiation that is centered 
about the tangent to the orbit at the source point. The observation point P 
and the vector r and n are therefore within the (y, z)-plane and radiation is 
emitted at angle 9 with respect to the z-axis. Following the above discussion 
the azimuthal angle is constant and set to <fi = ir/2. The definition of n and 
R can be used to compute the scalar product in the exponential term in (29) 
so that 



n 



R = R cos 9 sin ivot' , 



where 9 is the vertical angle and uoq = c/R. Since the angles are small and the 
relativistic 7 factor is large, we can replace 

cos 9 ~ 1 - 9 2 /2 , sin u t' ~ cu t' - ulit' f/Q . 



Therefore the exponential factor becomes 



R 

c 



f 



2/ + /\3 



f'--ccs^in^f'-^(l + 7 2 ^ 2 )+ ( . /; 



The triple vector product in (27) can be evaluated in a similar way as 

n x [n x 0} = u t'e x + 9e y , 

We can now write E r {u) as 



E T (u) = ^r^- < r, 



c r 



f 



uj t'exp iiu — (l+ 7 2 fl 2 ) + 



+e. 



00 / / / 

-rv-i * * 





(30) 



The integrals in (30) can be expressed in terms of modified Bessel functions 
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as first pointed out by Schwinger [4]. For this we need a change of variables 
as follows: 



r = 



(l + 7 202)l/2 



3/2 



where u c = 37 3 u> /2. With these substitutions we get integrals of the form 
°° 3 1 1 

/"-G«(-+r 3 )) dT= 73^ K) ' 



00 

/™(§< r+ r*)) dT = >■/■«>• 



The Fourier transform of the electrical field (30) finally becomes 



c r oj c 



, 7^i/3(0 
^2/3(0^ + i (1+72 g 2)1/2 e y 



Using (26) we get an expression for the spectral distribution: 

d 2 W 3e 2 iu 2 



7 2 (l+7 2 # 



2/)2\2 



^3(0 + 7^^/3(0 



l + 7 2 # 



(31) 



Quite often an observer is interested in the synchrotron radiation energy emit- 
ted over all vertical angles. Integration of equation (31), over the angle 9 gives 
the required result. We note that the solid angle d^ = d#d0 and that radi- 
ation does not depend upon the azimuthal angle. So we can write 



d 2 W 3e 2 u>V 



d^d< 



An 2 cuil 



J {i + i 2 e 



2/)2\2 



K 2/3(0 + I + ^202 



d9 



(32) 



The angle 9 appears in (32) in rather a complicated way which makes it diffi- 
cult to perform the integration directly. We shall not describe the integration 
process in detail, not only because it is available in the form of comprehensive 
treatises [3], but also because this is a pure mathematical problem. If we look 
it up in [3] we see that (32) can be written as 



d 2 W V3e 2 j 
dujd4> 47T 



OJc. 



K 5 /3(y)dy 



u>/u) c 



(33) 
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synchrotron radiation production from 
circular motion 



time variation of electric field energy spectrum 

at detector in far zone 




Fig. 13. Radiation field from an electron moving along a circle in the time and in 
the frequency domain 

radiation from an electron moving along the arc of a circle 

( L m > R/7 ) 



time variation of electric field energy spectrum 

at detector in far zone 
( electric field averaged over time / ) 




Fig. 14. Radiation field from an electron moving along an arc of a circle in the time 
and in the frequency domain 

Equation (33) is the required expression for the frequency spectrum of the 
radiation from an electron moving along a trajectory which is a circle. For a 
small argument w/w c C 1 we may apply an asymptotic approximation for the 
modified Bessel's function and get 

— / X 5/3 (y)dy^l.33 — for cu/u c < 1 . 

UJ C J W c / 

Ul/u> c 



Most textbooks on synchrotron radiation discuss equation (33). However, no 
attention is usually paid to the region of applicability of this equation. Calcu- 
lations which led to (33) are only valid for an electron moving in a circle. On 
the other hand we want to study synchrotron radiation from a dipole mag- 
net too. Consider what would happen if, instead of an electron moving in a 
circle, we had an electron which moves along the arc of a circle. Fig. 13-15 
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radiation from short bending magnet 

( L m « R/7 ) 



time variation of electric field 
at detector in far zone 
( electric field averaged over time / ) 



AtKw? = iV(2c7 ) 



energy spectrum 



Fourier 
transform 




Fig. 15. Radiation field from an electron moving along a short bending magnet in 
the time and in the frequency domain 

sketch the expected frequency spectrum. These figures provide an illustration 
of the transition from the usual synchrotron radiation by an electron in a 
circle to radiation by an electron moving along an arc of a circle. The first 
sketch (Fig. 13) describes radiation from a circle. It is useful to consider fur- 
ther the relationship between the time domain and the frequency domain. The 
frequency spectrum of the radiation pulse is given by the expression (33): the 
Fourier transform at uo = 0, E(0), is zero because E(t) averages to zero (see 
Fig. 3). 

In the case of dipole magnet, things are quite different. The frequency spec- 
trum may be obtained calculating the Fourier transform of the time distri- 
bution of the electric field at the observer's position. It can be obtained nu- 
merically without much work by noticing that we can use curves in Fig. 5-7. 
As a result, we can define general properties of the spectrum without any 
analytical calculations. Consider Fig. 14, 15. As these curves illustrate, when 
the electron moves along an arc of a circle, the spectral distribution does not 
tend to zero when uo — * 0. Although this fact may look surprising, this is quite 
natural, since the average of the electric field over time is nonzero when an 
electron moving along the arc of a circle. 



3.2.2 Radiation from an electron moving along the arc of a circle 

Many texts on synchrotron radiation theory do not seem to provide a deriva- 
tion of the expression for the radiation spectrum from an electron moving 
along an arc of a circle. Here we present such a derivation, which we believe 
to be quite simple and instructive. Analytical studies of this problem were 
first performed in [5]. A similar treatment may also be found in the book [9]. 
Our consideration is restricted to the simplest case from an analytical point 
of view, namely the case of a bending angle m very small compared with 
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Fig. 16. Geometry for synchrotron radiation from an electron moving along an arc 
of a circle 

the I/7. No assumptions on the magnetic field parameters have been made 
to derive the radiation spectrum in the form of equation (27): therefore, our 
starting point is the same as that for the radiation spectrum calculation from 
a circular motion in the previous subsection. Introducing the notation 

n x [{n — f3) x f3] 
(l-n-(3) 2 ' 



G(t') = 



we rewrite (27) in the form: 

00 

E T (u) = I G{t') exp[icu(t' - n ■ R(t')/c)]dt' , 

c 7"n J 



—00 



Now we have to perform the integration. The function G(t') is a smooth curve 
and does not vary very much across the very narrow bending angle region in 
the case of a short magnet: therefore we may replace it by a constant. In other 
words, we simply take G(t') outside the integral sign and call it G . In this 
case the integral over t' in the latter equation is calculated analytically 

- L m /(2c/3) 

E T (cu) = / exp[icu(t' - nR(t')/c)] df 

c r J 

-L m /(2c/3) 



(-e)G ^m sin 


[a- 


n0)ujL m /(2c(3) 


c 1 ft o 1 


a 


— ft 


■P)uL m /(2cP) 



(34) 
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Next we compute the quantities required for expressing explicitly the latter 
equation. If the system has no azimuthal symmetry, the angular distribution 
of radiation depends on two variables. These may be taken (Fig. 16) to be the 
angle 9 between n and the x — z plane, and the azimuthal angle 0. We can 

now write down the components of the vectors (3, j3, and n as: 

■ (3 2 c 
(3 = (3e z , (3 = (3e x = —e x , 

n = cos 9 sin (fie x + sin 9e y + cos 9 cos <pe z . 



These definitions can be used to get the triple vector product 



nx [(n- (3) x (3] = $ (cos 2 9 sin 2 - 1 + (3 cos 9 cos </>) e x 
+ (sin 9 cos 9 sin 0) e y + (cos 2 9 cos sin — (3 cos 9 cos 0) e z 



(35) 



We can use these expressions to compute the Fourier transform (34). The 
spectral distribution is proportional to the square of electric field and if we 
follow through the algebra we find that 

d 2 W e 2 L 2 m (Al + A\) sin 2 [(1 - (3 cos 9 cos <f>)uLj(2c/3)] 
dcudfi 4vr 2 i? 2 c(l -/3cos#cos0) 4 [(1 - /3cos#cos0)u;L m /(2c / 3)] 2 ' 



where the following notation has been introduced: A a = sin 0-/3 sin 9 and 
At, = cos 9 cos 0. The quantities A a , A n are characteristics of the a— and n— 
components of the single-particle linearly polarized radiation. 

We emphasize the following features of this result. The radiation is very 
much collimated in the forward direction, most of the energy being emitted 
within the cone AQ = A0A# of 7 -2 . The radiation maximum is reached 
at 9 — 0, = 0. Spectral properties of the radiation are defined by the 
function [sin 2 u/u c ]/[uj/uj c ] 2 , where the critical frequency is oj c = 2c/[(l — 
/3cos#cos0)L m ]. The spectrum from an electron in a "short" magnet is prac- 
tically "white noise" spreading from zero up to the frequencies oo ~ uo c . As one 
can see, the spectral properties of radiation considered here are significantly 
different from those of synchrotron radiation from an electron moving in a 
circle 3 . 



3 There is an obvious analogy between the synchrotron radiation from "short" 
magnet and radiation from bremsstrahlung effect. The mathematics of the two 
problems is essentially the same 
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3.3 Limitations of standard results 



As already pointed out, we can decompose the CSR spectrum, P(u), into 
the product of the square modulus of bunch form-factor | F(cu) | 2 with the 
radiation spectrum from one electron p(u). To find | F(u) | 2 from P(ou), one 
needs the quantity p(ou), which is usually very difficult to know with great 
accuracy within the long wavelength range. Actually, as already discussed, 
the well-known theory of synchrotron radiation uses approximations which 
are valid under certain conditions, inappropriate under others. In particular, 
Schwinger approach relies on several assumptions, which we review once again: 
first, the observer is located in such a way that the velocity term in the Lienard- 
Wiechert formula can be completely neglected and that the unit vector n can 
be considered constant throughout the electron evolution. Second, a circular 
trajectory is postulated. Finally, no aperture limitations is considered at all. 
These assumptions must be analyzed in order to understand how the CSR 
pulse is altered in realistic situation. In this subsection we will deal separately 
with all of them. 

Let us imagine that our electron bunch moves along an arc of a circle and 
that there is no aperture limitation. We can take into account, then, effects 
as a finite distance between the source and observer, a finite magnet length 
and the presence of velocity fields by means of (6) where the expression for 
the electric field is given by the strict Lienard-Wiechert formula. Analytical 
methods are of limited use in the study of CSR in the near zone, and numerical 
calculations must be selected. The application of similarity techniques allows 
one to present numerical results in such a way that they are both general and 
directly applicable to the calculation of specific device situations. The use of 
similarity techniques enables one not only to reduce the number of parameters 
but also to reformulate the problem in terms of variables possessing a clear 
physical interpretation. Each physical factor influencing the CSR production 
is matched by its own reduced parameter. For the effect under study this 
reduced parameter is a measure of the corresponding physical effect. When 
some effect becomes less important for coherent radiation, it simply falls out 
of the number of the parameters of the problem. 

For our present purposes we would like to concentrate completely on the 
temporal structure of the CSR pulse. For any CSR pulse there exist some 
characteristic time, which determines the CSR pulse duration. Its magnitude is 
of order of the inverse of the frequency spread in the CSR beam. The behavior 
of the CSR pulse profile as a function of dimensionless parameters provides 
information on the spectrum distortion. To be specific, we consider the case of 
an electron beam with a Gaussian axial profile of the current density. When 
comparing the temporal structure of CSR pulses, it is convenient to use a 
normalized field amplitude: here, the normalization is performed with respect 
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to the maximal field amplitude. The rms electron pulse duration is assumed 
to be large <7t ^> R/(c~y 3 ) and we focus on the radiation pulse seen about the 
tangent to the orbit at the middle point of magnet. When the vertical angle 
9 = the normalized coherent field amplitude E x (t) / E max is a function of six 
dimensional parameters: 

t , v , R , L m , <7 T , | r | . 

It is relevant to note that only two dimensions (length and time) are sufficient 
for a full description of the field profile: after appropriate normalization it is 
a function of four dimensionless parameters only: 

E x {t)/E max = E X = D(t, a, p, f ) , (36) 

where t — tju^ is the dimensionless time, cr = c^oO"t7 3 is the dimensionless 
electron pulse duration, | f | (ca^^^R^ 2 ^ is the dimensionless distance 
between source and observer, p = 0^/ (Gc^o^t) is the magnet length parameter. 
In the general case the universal function D should be calculated numerically 
by means of strict Lienard-Wiechert formulas. 

The changes of scale performed during the normalization process, mean that 
we are measuring time, bunch length, distance from the source and magnet 
length as multiples of "natural" CSR units. There is a physical reason for being 
able to write the field profile as in (36): let us explain this fact beginning with a 
qualitative analysis of the radiation from an electron moving in a circle, in the 
long wavelength asymptotic. Synchrotron radiation is emitted from a rather 
small area and we need to determine this area for observers whose detection 
systems collect information over a long time period or- The radiation pulse 
length is equal to the time taken for the electron to travel along any arc AB, 
reduced by the time taken for the radiation to travel directly from the arc 
extreme A to B. Between point A and point B we have a deflection angle 
0, so that o"t — R(f)/((3c) — 2_Rsin(0/2)/c, and sin(0/2) can be approximated 
by sin(0/2) ~ 0/2 — 3 /48 for small angles. Then the pulse duration reduces 
to a T ~ i?0 3 /(24c). The radiation source extends over some finite length 
R(f) ~ L { — (ca T y/ 3 R 2 / 3 along the particle path. We see that the reduced 
distance can be expressed as fo =| r*o | /Lf. One can find that the ratio 
(i?0 m ) 3 / (6L 3 ) is equal to p, which we use now as a measure of finite magnet 
length effects. 

3.3.1 Finite distance effects 

Let us imagine that our electron bunch moves in a real circle and there is 
no aperture limitation. We can first consider the contribution given by the 
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Fig. 17. Time structure of a CSR pulse from a Gaussian electron bunch moving 
in a circle at different reduced distance between source and observer. Here 6 = 0. 
The dashed curve is calculated within the far zone approximation (19). Numerical 
calculations have been performed with strict Lienard-Wiechert formula 

acceleration field alone and then focus on the contributions by the velocity 
field. When the electron bunch radiates from a circle the spectrum is obviously 
independent of the p parameter. In the long wavelength asymptotic (a ^> 
1) the acceleration field is described by one dimensionless parameter only: 
f , where f =| f | (cot) -1 ' 3 ^ -2 ' 3 is the dimensionless distance between 
radiation source and observer. 

The region of applicability of Schwinger's formulas requires the dimension- 
less distance to have a large value f ^> 1. In fact, as previously explained, 
the radiation source extends over some finite length Rip ~ Lf, and this length 
corresponds to a transverse size of the radiation source d ~ Rcj) 2 . The vector n 
changes its orientation between point A and point B of an angle of order d/r , 
where r$ is the distance between source and observer. Our estimates show 
that in the case d/ro <C <p, the vector n in Lienard-Wiechert formula is almost 
constant when the electron moves along the formation length Rip. Thus, we 
can conclude that the unit vector n can be considered constant throughout 
the electron evolution only if r Lf. 

The results of numerical calculations for several values of f o are presented in 
Fig. 17. Calculations have been performed using the strict Lienard-Wiechert 
formula. The plots in Fig. 17 give an idea of the region of validity of the far 
zone approximation considered above. At large distance the CSR pulse profile 
is simply the far zone profile (19). One can see that (19) works well at f = 100. 



32 




Fig. 18. Illustration of similarity techniques. Time structure of a CSR pulse from 
a Gaussian electron bunch moving in a circle for various sets of parameters. The 
reduced distance is held constant, fo = 15. The solid curve is calculated within 
the far zone approximation (19). Numerical calculations have been performed with 
strict Lienard-Wiechert formula 
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Fig. 19. Illustration of similarity techniques. Time structure of a CSR pulse from 
a Gaussian electron bunch moving in a circle for various sets of parameters. The 
reduced distance is held constant, fo = 0.15. Numerical calculations have been 
performed with strict Lienard-Wiechert formula 

Then, at f = 3, the CSR pulse envelope is visibly modified. As the distance is 
decreased, the difference between the approximate and the strict pulse profile 
becomes significant. From a practical point of view this set of plots covers all 
the region of interest for the distance between observer and sources. 
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To check that no mistakes have been made in our similarity techniques we 
evaluate the normalized CSR pulse profile for several sets of problem param- 
eters. The reduced distance is held constant. In the present example, these 
values are r = 0.15, 15. The plots are calculated, numerically, using the strict 
Lienard-Wiechert formula. Fig. 18, Fig. 19 show such plots. It is seen that 
there is a good agreement with the prediction that the acceleration field from 
a circle in long wavelength asymptotic is a function of the reduced distance 
only. 

3.3.2 Velocity field effects 

Usual theory of synchrotron radiation is based on the assumption that ac- 
celeration field dominates and all the results presented above refer to the case 
when there is no influence of the velocity field on the detector. The acceler- 
ation field dominates in the far zone only, and we want to study near zone 
effects too. The physics of coherent effects studied by means of general expres- 
sions for the Lienard-Wiechert fields is much richer than that of the simplified 
radiation field model considered above. In the long wavelength asymptotic the 
velocity part of the coherent electric field from a particle in a circle is a func- 
tion of two dimensionless parameters: the reduced distance parameter f and 
the reduced electron pulse duration parameter a: 

£(v) (t)/£max = E {v) = D(t, a, f ) , (37) 

where the normalization of the velocity field is performed with respect to the 
maximal acceleration field amplitude. 

To show that this is correct we can perform a simple check. In Fig. 20 we 
present numerical calculation results for the velocity field for different sets of 
parameters. The reduced distance and the reduced electron pulse duration are 
held constant. It is seen that all numerical results agree rather well. Therefore 
the plots in Fig. 20 convince us that the dimensionless equation (37) provides 
an adequate description of the numerical calculations. 

Fig. 21 and Fig. 22 show the dependence of the normalized velocity field 
amplitude on the value of the reduced electron pulse duration for different 
values of the reduced distance. Using these plots we can give a quantitative 
description of the region of applicability of the radiation field model. It is seen 
from the plots in Fig. 22 that in the near zone we cannot neglect the influence 
of the velocity field on the detector. 

Let us now estimate the importance of the velocity field effect. Let us start 
considering the total velocity field pulse as a superposition of single particle 
fields at a given position in the far zone. To calculate the integral one should 
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Fig. 20. Illustration of similarity techniques. Electric field pulse due to velocity 
term from a Gaussian electron bunch moving in a circle for various sets of param- 
eters. The reduced distance and reduced bunch length are held constant, fo = 1 
and ctt = 100 respectively. Numerical calculations have been performed with strict 
Lienard-Wiechert formula 

take into account the property of its kernel. In the far zone the velocity field 
from one electron is close to an antisymmetric function and the average of 
the electric field over time is close to zero. The approach used in section 
3.1 can be also used to compute the coherent velocity field. Under the long 
electron bunch condition the kernel (velocity field from one electron) can be 
substituted by its asymptotic behavior. If we wish to estimate the normalized 
amplitude of the coherent velocity field we can get it simply by dividing the 
asymptotic behavior of the velocity field kernel by the asymptotic behavior of 
the acceleration field kernel, so that, in the far zone, the normalized velocity 
field is of order E v /E acc ~ R/('y 2 9 2 r ), where 9 ~ (ca^/R) 1 ^ is the natural 
synchrotron radiation opening angle with frequency uj ~ a^ 1 <C cy 3 /R. Using 
normalized variables we get 

E (v) ~ r^fr- 2 ^ for <t, f > 1 . 



As we can see from Fig. 21, numerical calculations in the far zone confirm 
this simple physical consideration. The value of E v that we can expect is 
found remembering that, in the example given in Fig. 21, f = 100, cr = 100; 
therefore the normalized velocity field would be about 0.0004, which is the 
same order of magnitude as results of numerical calculations (0.0002). Also, 
E v varies roughly as a~ 2 / 3 . 
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Fig. 21. Electric field pulse due to velocity term from a Gaussian electron bunch 
moving in a circle at different reduced bunch length. The reduced distance is held 
constant, rg = 100. Here 6 = and the normalization is performed with respect to 
the maximal acceleration field amplitude 

The normalized velocity field amplitude decreases, as we see from our esti- 
mations, linearly with distance, which means that if we are in the near zone 
at r = 1 there will be E v ~ 0.04. Numerical calculation results show that, 
at the value f = 1, we have E v = 0.2. We should say that our approximate 
treatment of coherent velocity field breaks down once source and observer get 
as close as they are at r ~ 1. We have, however, a simple explanation for that. 
In addition to the antisymmetric part of kernel we have just described, there 
is also a symmetric part. When source and observer are far apart the observer 
sees only the antisymmetric part and the average of electric field from one 
electron is close to zero. At very close distances there begins to be some extra 
symmetric part of the kernel. This symmetric field, which also varies with the 
separation, should be, of course, included in more precise estimates. 



3.3.3 Finite magnet length effects 

Until now we have considered the case in which electrons moves in a circle. 
Here we will generalize this situation and study the case where the electrons 
move along an arc of a circle. In order to fully characterize the magnet length 
required to assure the accuracy of the circular motion model, it is necessary 
to specify the distortion of Schwinger spectrum associated with all interesting 
magnet lengths in an experiment. Normalizing equation (24), we obtain that, 
under the far field approximation, the CSR pulse profile is a function of only 
one dimensionless parameter p = 0^/(6u;o0t), where m is the magnet angu- 
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Fig. 22. Electric field pulse due to velocity term from a Gaussian electron bunch 
moving in a circle at different reduced bunch length. The reduced distance is held 
constant, fo = 1. Here 6 = and the normalization is performed with respect to 
the maximal acceleration field amplitude 

lar extension. The region of applicability of Schwinger's formulas require the 
parameter p to have a large value p 1. Using the plots presented in Fig. 23, 
one can characterize quantitatively the region of applicability of the circular 
motion model. 



3.3.4 Diffraction effects 

Another problem concerns radiation spectrum distortions: this has to do 
with aperture limitations. In real accelerators, the long wavelength synchrotron 
radiation from bending magnets integrates over many different vacuum cham- 
ber pieces with widely varying aperture. The SR spectrum formed by such a 
system will be perturbed by the presence of an aperture and we should seek 
a mathematical mean for predicting these perturbations. 

One may wonder how the existence of an aperture in an accelerator vacuum 
chamber and in a beam line, which is relatively far from the source, can influ- 
ence the synchrotron radiation spectrum. To illustrate this point, we consider 
the situation shown in Fig. 24. The object of interest is an electron moving 
in a circle. Between the observer and the source there is an aperture with 
a characteristic dimension d. Qualitatively, an observer looking at a single 
electron is presented with a cone of radiation characterized by an aperture 
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Fig. 23. Time structure of a CSR pulse in the far zone from a Gaussian electron 
bunch moving along an arc of a circle at different reduced magnet lengths. Here 
6 = 0. Calculations have been performed with anaytical formula (24). Results at 
p = 50 (solid curve) agree with our analytical expression (19) for circular motion 
within 1% 



angle of order 9 ~ yd/R. Fig. 24 shows part of the trajectory of an electron 
travelling along an arc of a circle of radius R. At first glance this intuitive 
argument indicates that the frequency spectrum may be obtained by reducing 
the situation shown in Fig. 24 to Fig. 4 and simply calculating radiation from 
an electron moving along a trajectory which is an arc of a circle. However, the 
situation is more complicated than that. The field at P can be represented as 
a sum of two parts - the field due to the source plus the field due to the wall, 
i.e. due to the motions of the charges in the wall. The latter kind of radiation 
field can be described in terms of diffraction theory. Diffraction is the process 
by which radiation is redirected near sharp edges. As it propagates away from 
these sharp edges, it interferes with nearby non-diffracted radiation, produc- 
ing interference patterns. A finite aperture introduces, therefore, diffraction 
effects specific to the geometry and clearly dependent on the wavelength. For 
structures such as pinholes it is found that these diffraction patterns propagate 
away from the structure at angles of order 9& — c/ ud, where d is the character- 
istic aperture dimension. The region of applicability for the far diffraction zone 
where effects play a significant role is given by the relation L a ^> L d ~ ud 2 /c, 
where L a is the distance between observer and aperture. When the wavelength 
is about R9 3 ~ (rf 3 //?) 1 / 2 , the latter condition transforms to L a ^> L d ~ R9. 
This simple physical argument indicates that the situation shown in Fig. 24 
at L a > R9 ~ (Rd) 1 / 2 cannot be reduced to the situation shown in Fig. 4. The 
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Fig. 24. Depiction of the effects of an aperture limitation. CSR from a circle allows 
us to measure the bunch profile form-factor in a frequency range down to cj a only- 
significance of the discussed limitation cannot be fully appreciated until we 
determine typical values of the parameters that can be expected in practice. 
For example, if R = 3m and A = 100/im, 9 ~ 0.03, so that at distance greater 
than 10 cm from the aperture, a 1 cm diameter hole would be sufficiently 
narrow to perturb significantly the synchrotron radiation spectrum at that 
wavelength. 

The computation of spectrum distortions for practical accelerator applica- 
tions is a rather difficult problem. One encounters vacuum chamber compo- 
nents which require special attention and sophisticated techniques should be 
developed to deal with this situation. The character of spectrum distortions 
largely depends on the actual geometry and on the material of the vacuum 
chamber, so that we may expect significant complications in the determination 
of these distortions along a vacuum enclosure of the actual device. The vacuum 
chamber of an accelerator presents a too complicated geometry to allow an 
analytical expression for spectrum distortions. In principle each section of the 
vacuum chamber must be treated individually. By employing two or three- 
dimensional numerical codes it may be possible to determine the spectrum 
for a particular component. Yet, every accelerator is somewhat different from 
another and will have its own particular spectrum distortions. For this reason, 
in this discussions we focused only on basic ideas about aperture limitations. 



4 Radiation from an undulator 

To optimally meet the needs for beam diagnostic measurements with the 
features of synchrotron radiation, it is desirable to provide specific radiation 
characteristics that cannot be obtained from bending magnets but require 
special magnetic systems. To generate specific synchrotron radiation charac- 
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teristics, radiation is often produced by special insertion devices installed along 
the particle beam path. An insertion device does not introduce a net deflec- 
tion of the beam and we may therefore choose any arbitrary field strength 
which is technically feasible to adjust the radiation spectrum to experimental 
needs. In this Section we discuss the basic theory of one such devices, a planar 
undulator, and its applications to CSR diagnostics. 

4-1 Basic theory of undulator radiation 

The magnetic field on the axis of a planar undulator is given by 

H(z) = e y H w sin(k w z) , 

where e y is the unit vector directed along the y axis of the Cartesian coordinate 
system (x, y, z) in Fig. 25. The Lorentz force is used to derive the equation of 
motion of an electron with energy / -fm c c 2 in the presence of a magnetic field. 
Integration of this equation gives 

v x (z) = e x c9 w cos(k w z) , 

where 9 W = if/7 and if = eH w / (m e c 2 k w ) is the undulator parameter. It is 
useful to present another form of this expression convenient for numerical cal- 
culations: if = 0.093 x ii w (kGs) x A w (cm) where A w = 2n/k w . The transverse 
and longitudinal velocity components are related by (z) 2 = j3 2 c 2 — (x) 2 , which 
gives 

K 2 

Pz = Aw - cos 2k w z , 

where Aw = (3[l — K 2 /(A^ 2 )] is the mean velocity of the electron along the 
undulator. 

It is generally assumed that the distance between the undulator and the 
observation point is much larger than the undulator length. This assumption 
is usually referred to as far field approximation. Suppose this condition is met. 
Denoting with r, as usual, the vector from the electron position at the retarded 
time to the observation point and with r the vector from the coordinate origin, 
placed somewhere inside the undulator, to the observation point we compute 
r ~ f, 3 — n ■ R with n = r / \ r \. The unit vector from the retarded source 
particle to the observer is, from Fig. 25 

n = cos sin 9e x + sin <f> sin 6e y + cos 9e z . 
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Fig. 25. Radiation geometry for the undulator field 
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Fig. 26. Synchrotron radiation production from successive poles 

The vector R, which describes the motion of the electron in the x — z plane, 
can be obtained from the previous treatment: 



K K 2 
x(t) = - — sin(k w p av ct) , z(t) = (3 av ct - — 



■ sin(2/c w/ 9 av ct) . 



Let us find the Fourier harmonics of the radiation field from an electron moving 
along the undulator. As in section 3.2, our starting point is equation (27) since, 
once again, this is a rigorous equation under the assumption of a large distance 
between observer and radiating electron. Note that the integration in (27) is 
to be performed over all time but in practice, without loss of generality, it 
can be restricted to the time during which the electron is passing through the 
undulator. We can write: 



HO 

c | r o 



nx[(n — P) x (3} 
(1 -n-p) 2 



exp[i cu(t' - n ■ R{t')/c)\ dt' ,(38) 
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where T = A w / (c/3 av ), and iV w is the number of undulator periods. The integral 
in the latter expression can be represented as a sum of iV w integrals over 
each undulator period (see Fig. 26). Integrating by parts (38) and taking into 
account the sum of the geometric progression 

£ exp(in0) = ""y/ 2 exp i(7V w + 1)0/2 



n=l 



sin 0/2 



we find 



-± , . (— e) n x In x 3] 
E r (u) = -^/j { ^ ^ exp 



c\r \ (l-n-0) 

s inAU>(l-/3 av n 2 ))T/2 
sinu;(l - 3 w n z )T j2 



\uj{t' -n- R{t')/c) 



t'=N w T 



t'=0 



iw- 



exp[i(iV w -lV(l-/5 av n 2 )T/2] 



x Jnx[nxf3}exp \uj{t' — n • R(t')/c) dt' 



(39) 



For our purposes, it is most convenient to express the latter equation in the 
form 



2i(-e) 



c | f | (1 -/3 av cos6>) 
1 - /3 av cos 6>) 



exp [i(iV w - l)w(l - /? av cos 0)T/2] 



x { n x [n x 8 ] 

I (1-n-A)) 

x exp [iw(l - /3 av cos#)T/2] sin [JV w a;(l - /3 av cos 0)T/2] 

smN w u{l- 3 &v cos6)T/2 



-[w(l-0 w cos9)T/2] 



sinu;(l — 3 W cos9)T/2 



T 



x— nx [nx f3]exp\iuj(t' -n- R(t')/c)}dt'\ . 



(40) 



Here cfio = c(3(0) = c3(N w T) is the initial and final velocity of the elec- 
tron. From the shape of the interference function sin 2 (N w x)/ sin 2 (a;) we pre- 
dict that important contributions to the second term will occur when x = 
uj(1 — 3^-v cos 0)T/2 does not differ significantly from the resonance value x = ir. 
The maximum of the absolute value at resonance is, then, proportional to the 
number of periods N w and, on the assumption of N w ^> 1, we can neglect the 
first term. 
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The spectral distribution of the radiation energy is given by equation (26). 
Taking the square modulus of (40) and keeping the resonance term only we 
find 



sin iV w c<;(l - f3 &v cos 9)T/2 
sinu;(l — f3 av cos9)T/2 



d 2 W _ eV 
d 00 d 47r 2 c 

T 

x /nx[nx/3]exp ioo(t' - n ■ R(t')/c) dt' 



(41) 



Let us now analyze the interference function in (41) 



siniV w w(l - (3 av cos 9)T/ 2 
sincj(l - (3 av cos9)T/2 



and study some of its consequences. The whole curve is localized near 

-l 



u; = c<j rcs = 2ck w Y 



K 2 
2 ' 



The latter condition is called resonance condition and tells us the radiation 
frequency as a function of the undulator period A w , the undulator parameter 
K w , the electron energy / ym e c 2 , and the polar angle of observation 9. Let 
us concentrate on the radiation produced in the forward direction. In this 
case 9 = and there is no vertically polarized radiation. By means of an 
angular filter, whose principle is shown in Fig. 27, we can obtain the undulator 
radiation from the central cone. Note that for the radiation within the cone 
of half angle 



^cen 



Kl/2 



the relative spectral FWHM bandwidth is Au/w = 0.89/iV w . This means that 
we have a very sharp central maximum with very weak subsidiary maxima 
on the sides. In fact, it is possible to prove that the intensity at the next 
maximum is less than 5% of the first maximum intensity (see Fig. 28). There 
is another important feature of the interference function: if the frequency is 
increased by any multiple of w res , we get other strong maxima of interference 
function at uj = 2u res , 3u res , 4u; rcs , and so forth. Near each of these maxima 
the pattern of Fig. 28 is repeated. 
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Fig. 27. Undulator radiation with an angular filter 
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Fig. 28. Radiated wavetrain from a single electron traversing an undulator, as de- 
tected into the central cone, and corresponding spectral distribution function 

There is no forward radiation for even harmonics of the fundamental fre- 
quency u res . Let us briefly justify the latter statement. First let us calculate 
the integral in (41). At 9 = the triple vector product is 



n x [n x (3\ = (3 — n(n ■ (3) = f3 x e x = — cos(/c w z)e a 

7 



(42) 



The exponential factor in (41) can be evaluated in a similar way 



cxp 



i muJo(t' — ft • R{t') / c) = exp [imk w z + iQsm2k w z] , 



(43) 



where m is integer, 

Q = K 2 muJo/(8k wl 2 ) = mKl/(4 + 2K 2 W ) , 
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and u = 2 7 2 A; w /[c(l + Kl/2)} is the resonance frequency of the fundamental 
harmonic at 6 — 0. Using (42) and (43), together with the expansion 

n=+oo 

exp[iQsin(2A; w 2;)] = ^ J n (Q) exp(2 ink w z) , 



where J n is the Bessel function of order n, we obtain: 
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n= J rOO 



X! Jn(Q) / cos(k w z) exp[i(2n + m)k w z] d- 



(44) 



which is non-zero only when m + 2n ± 1 = 0. When m is odd the spectral 
and angular density of the radiation energy emitted at zero angle by a single 
electron during the undulator pass is given by the expression 



d 2 W(uj m ) 



00,, 



6»=0 



c(l + K 2 /2) 2 [^(a;-^)/^ 



where the following notation has been introduced: uj m = mooo, 

A,U = [J(m-l)/2(<5) - J(m+1)/2(Q)] ■ 



When deriving the latter expression, we used the following relation for the 
Bessel functions: J- n (Q) = (—l) n J n (Q). Now we can proceed to calculate 
the energy radiated into the central cone at the fundamental harmonic. In 
the small angle approximation the solid angle is equal to d^ = # d # d </>• 
Integration of spectral and angular density over u and <fi gives us factors ujq/N w 
and 2n respectively. We also have integrate over 9 from to 9 cen . Altogether, 
the energy radiated into the central cone by a single electron is given by 

CCD " c(l + i^/2) ' 145 ] 



where, in this case, A u = [J (Q) - J^Q)} and Q = fT 2 /(4 + 2Kl). 
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Fig. 29. Scheme for electron bunch length diagnostics based on CSR from an undu- 
lator 

4-2 Bunch length diagnostic technique based on coherent radiation from an 
undulator 



Above, a brief treatment of undulator radiation was given in view of beam 
diagnostics applications which we are going to describe in this paragraph. Our 
technique relies on using an electromagnetic undulator and on recording the 
energy of the coherent radiation pulse into the central cone. The configuration 
is sketched in Fig. 29. The electron bunch passes through the undulator and 
produces a CSR pulse at some specific resonant frequency ujq{K). The signal 
(energy of the CSR pulse within central cone) is recorded by a bolometer; 
as we will explain in this subsection, repetitions of this measurement with 
different undulator resonant frequency allow one to reconstruct the modulus 
of the bunch form factor. 



The undulator has, in our scheme, a high magnetic field and a large K >> 1 
value. A 40 cm magnetic period and a relativistic factor 7 ~ 10 3 would produce 
a fundamental wavelength of about 100 microns when K ~ 30. On the other 
hand, the transverse velocity of the electron is small: j3 x <C 1. The combination 
of these conditions gives an inequality 1 <C K <C 7 which, with a little loss of 
generality, we will consider satisfied. 

The principle of operation of the proposed scheme is essentially based on 
the spectral properties of the undulator radiation. A sample of undulator 
radiation spectrum at zero angle is shown in Fig. 30. The distribution of the 
radiated energy within different harmonics depends on the value of K. For 
A" > 1 we will get strong undulator maxima at 00 = 3ujq, 5uq and so forth (as 
already discussed, the radiation spectrum along the undulator axis includes 
only odd harmonics). The bunch form factor is the exponential function (1) 
and the square modulus of bunch form-factor falls off rapidly for wavelengths 
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Fig. 30. Some typical spectra of the longitudinal electron distribution and un- 
dulator spectral distribution function. The bandwidth of the undulator radiation 
(Alu/loo = 0.89/iV w ) must always be as small as the desired spectral resolution 



shorter than the effective bunch length. The fast drop off is evident and for 
wavelength of about three times shorter than the effective bunch length the 
radiation power is reduced to about 1 % of the maximum pulse energy at the 
fundamental harmonic. As a result, sharp changes of the bunch form-factor 
result in an attenuation of the high undulator harmonics. This point is made 
clear upon examination of the spectra presented in Fig. 30: actually, the bunch 
form-factor cut-off affects all higher undulator harmonics which are, therefore, 
suppressed. 

Our goal is to find the relationship between the expected value of radiation 
energy per pulse and the square modulus of the bunch form-factor. In general, 
| F(uj) | 2 will vary much more slowly in oo than the sharp resonance term. The 
two functions might appear as shown in Fig. 30. In such cases, we can replace 
| F(uj) | 2 by the constant value | F{uj q ) | 2 at the center of the sharp resonance 
curve and take it outside of the integral. What remains is just the integral 
under the curve of Fig. 28, which is, as we have seen, just equal to (45). Thus, 
the energy radiated into the central cone by an electron bunch is given by 



We emphasize the following feature of this result: the CSR radiation energy 
per pulse is proportional to the square modulus of the bunch form-factor at 
the resonant frequency. This fact allows one to reconstruct such quantity by 
repeated measurements of APFcsr at different resonant frequencies ojq. 



A^csr 




N 2 | F(u ) 



2 



(46) 
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4-3 Physical remarks 

The typical textbook treatment of the undulator radiation problem (see [3], 
[8]) consists in finding the expression for the undulator radiation spectrum 
from (29). However, no attention is usually paid to the region of applicability 
of this expression. In general this is not correct. In the undulator case, it is 
important to realize that (29) is valid only when the resonance approximation 
is met. Naturally, the resonance approximation cannot be expected to be valid 
when an insertion device has only a few periods. A second important limitation 
concerns the far field approximation. It is generally assumed that the distance 
between the undulator and the observer is much larger that the undulator 
length. This approximation is often assumed without any real justification 
other than the simplification that results. Such an assumption may or may 
not be valid, depending on the case. 



5 Constrained deconvolution 

If both modulus and phase of the bunch form-factor were somehow mea- 
sured, we would then known the Fourier spectrum of the bunch; an inverse 
Fourier transform of the measured data F(iv) would then yield the desired 
profile function F(t) with a resolution limited by the maximum achievable 
frequency range. Unfortunately, in practice, it is impossible to extract the 
phase information from the CSR measurement. A more realistic task is to 
measure the modulus of bunch form-factor only: in Section 4 we proposed a 
method to perform such a measurement, but independently on the technique 
selected, the following problem arises, to study what information about the 
bunch profile can be derived from the collected experimental data: the prob- 
lem is generally referred to as phase retrieval. There exist some cases in which 
the absence of phase information is of no consequence. For example, any sym- 
metrical density distribution can be recovered from the knowledge of the form 
factor modulus only. An example is given by a bunch with a Gaussian density 
distribution. In practice, the situations in which this assumption is satisfied 
are probably rather limited. 

A possible solution to the loss of phase information was suggested by Wolf 
[10]. A step forward toward the solution of the problem, is based on the fol- 
lowing consideration: if 

F(u) =| F(u) | exp[i0(cj)] , 
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then 



In 



F(u) = In | F(u) | + i</>(u) . 



It can be proved that if In 
recoverable from the modu 



F(u) is an analytical function, the phase will be 
us by the relation 



2lu ? ^ I F(u') 



Unfortunately, the analyticity of F{ui) is not a sufficient condition to assure 
analyticity for In F{uj) in that same region. The most obvious reason for 
this is the possible existence of zeros of F(uj) which lead to singularities in 
In F{uj) . In some situation the function F{uj) may have no zeros, in which 
case the so-called minimum phase solution is valid [11]. However, in general, 
zeros will be present, and their localization will be unknown a priori. 

It should be noticed that, whereas the dispersion relation method is very 
elegant and conceptually simple, it is not necessarily the "simplest" method to 
use in practice. For any given problem, the various possible approaches should 
be considered, for one may be distinctly easier than another, depending on 
the problem at hand. 

In the particular case of interest, namely an electron bunch originating by 
the XFEL driver accelerator, we can combine a-priori knowledge about the 
bunch density distribution with measured information about the bunch form 
factor. We know that, in the case under examination, the bunching process in 
zero order approximation can be treated by single particle dynamical theory. 
We cannot neglect the CSR induced energy spread and other wake fields when 
we calculate the transverse emittance, but we can neglect them when we cal- 
culate the longitudinal current distribution, simply because their contribution 
to the longitudinal dynamics is small. As a result, we will show that we can 
define a general temporal structure of the electron bunch after compression 
without any measurements. The purpose of the measurement is then to deter- 
mine the numerical value of the parameters on which such temporal structure 
depends. In this case, sufficient information is provided by the modulus of the 
bunch form-factor, allowing us to ignore the phase information. The process 
of finding the best estimate of the profile function F(t) for a particular mea- 
sured modulus of the bunch form-factor, including utilization of any a priori 
information available about F(t) refers to a technique called "constrained de- 
convolution" which will be exemplified below in detail. Although the context in 
which the suggestion of this method was made first was spectroscopy [12], the 
ideas apply equally well in the present context of bunch profile reconstruction. 



49 



410 




390^ . , . , . 1 \ , , , 1 , , , 1 , 1 

10 20 30 0.0 0.2 0.4 0.6 0.8 1.0 

t [ps] t [ps] 



Fig. 31. Typical phase space distribution of electrons after full compression with a 
single bunch compressor. The plot on the right presents an enlarged fraction of the 
plot on the left. The head of the bunch is at the left side 



Let us consider an example which shows the deconvolution process in easy- 
to-understand circumstances. For emphasis we choose, here, an experiment 
aimed at measuring the strongly non-Gaussian electron bunch profile at TTF, 
Phase 2 in femtosecond mode operation. The femtosecond mode operation is 
based on experience obtained during the operation of the TTF FEL, Phase 1 
[6] and it requires one bunch compressor only. An electron bunch with a sharp 
spike at the head is prepared, with rms width of about 20 microns and a peak 
current of about one kA. This spike in the bunch generates FEL pulses with 
duration below one hundred femtoseconds. An example of the longitudinal 
phase-space distribution for a compressed beam with RF curvature effect is 
shown in Fig. 31, when the longitudinal bunch charge distribution involves 
concentration of charges in a small fraction of the bunch length. In the fem- 
tosecond mode the first TTF magnetic chicane only will be operational, and 
this will be the default mode of operation for first lasing. 

Bunch compression for relativistic situations can effectively be achieved, at 
the time being, only by inducing a correlation between longitudinal position 
and energy offset with an RF system and taking advantage of the energy de- 
pendence of the path length in a magnetic bypass section (a so-called magnetic 
chicane) . In order to study the CSR diagnostic for a compressed bunch, let us 
first find the longitudinal charge distribution for our bunch model when it is 
fully compressed by a chicane (here we assume that the perturbation of the 
bunching process by CSR is negligible). Downstream of the chicane the charge 
distribution is strongly non-Gaussian with a narrow leading peak and a long 
tail, as it is seen from the simulation results in Fig. 31. Our attention is fo- 
cused on determining the longitudinal density distribution after compression. 
As first was shown in [13] there is an analytical function which provides quite 
a good approximation for the simulation data. The current distribution along 
the beam after a single compression is given by 
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Fig. 32. Current distribution along the bunch reconstructed from an ideal set of 
experimental data for the form-factor and correct value of time constant for the 
bunch tail (solid curve). Circles present data from a tracking simulation code 

I(t) ~ Joexp ( — —2 ] for ti > t > — oo ; 



2r 2 
M) 

t + t )/n 



I(t)c^tM for t> tl >0. (47) 



Equation (47) gives a good fit. To illustrate how good an approximation to 
the actual current distribution is provided by (47), we show in Fig. 32 the 
current as a function of time, first the numerical simulation data, then as 
computed from (47). Equation (47) involves four independent parameters in 
the calculation of longitudinal density distribution: r , which characterizes the 
local charge concentration (non-Gaussian feature) of the compressed density, 
ti, which characterizes the bunch tail and fitting parameters, to and t±. In 
this case detailed information on the longitudinal bunch distribution (actu- 
ally constants: ro, T\ and to, ti) can be derived from the particular measured 
square modulus of the bunch form-factor by using an independently obtained 
tail constant T\ (for example by means of a streak camera). A quantity of con- 
siderable physical interest is the parameter r which characterizes the leading 
peak. The value of r can be related to the value of the initial local energy 
spread A7/7 in the electron beam if the compaction factor i? 56 of the mag- 
netic chicane is known. Fig. 33 shows plots of the peak current versus time 
for several values of the local energy spread. By comparing these curves, some 
feeling can be obtained about the effective reduction in peak current. As the 
energy spread is increased from 5 keV to 15 keV, the maximum peak current 



51 




t[ps] 

Fig. 33. Typical current distribution along the bunch after full compression with a 
single bunch compressor. Curve 1, 2, and 3 correspond to local energy spread of 5, 
10, and 15 keV at the entrance of the bunch compressor. The head of the bunch is 
at the left side. The charge of the bunch is 3 nC 

falls from 3.5 kA to 2 kA. 

The highest time resolution obtainable by means of a streak camera is in the 
sub-picosecond range. By means of coherent radiation, the resolution can be 
made comparatively high because the resolution is determined by the spectral 
range in the measurement. Hence, this method will be useful for monitoring 
an ultra-short electron bunch at TTF FEL. In the TTF FEL Phase 2 the mea- 
surement of CSR spectrum can be made by using long period electromagnetic 
undulator as shown in Fig. 34. As already said, for the complete determina- 
tion of the bunch profile function, downstream of the first TTF FEL bunch 
compressor, a measurement of the bunch length can be made independently 
by using a sub-picosecond resolution streak camera. The limiting resolution 
of the streak camera will prevent the measurement of the spike width. Never- 
theless, the streak camera can detect the long time tail which is outside the 
CSR detector response. 

The significance of the proposed scheme cannot be fully appreciated until we 
determine typical values of the energy in the radiation pulse detected within 
the central cone that can be expected in practice. In the case of TTF FEL, 
the bunch form factor modulus falls off rapidly for wavelengths shorter than 
40 microns. At the opposite extreme, the dependence of the form factor on the 
exact shape of the electron bunch is rather weak and can be ignored outside 
the wavelength range 40/im < A < 100/im. This point is made clear upon 
examination of the curves presented in Fig. 35. The energy in the radiation 
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Fig. 34. Bunch length measurement scheme for TTF XFEL. The CSR source uses 
the spent electron beam coming from the XFEL 
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Fig. 35. Form- factor of the electron bunch corresponding to the bunch profiles shown 
in Fig. 33 

pulse can be estimated simply as in (46). From Fig. 36 it is quite clear that in 
the 40/iin — 100/xm wavelength region of the CSR spectrum any local energy 
spread in the electron beam smaller than 20 keV produces energies larger than 
0.1/xJ in the radiation pulse. 

To conclude this section we present examples of deconvolution of computer 
generated data; we hope that these will provide useful insight into the practical 
limits of deconvolution. We will demonstrate that constrained deconvolution 
is a useful experimental tool testable at each step in its development. By far 
the easiest and most informative way to examine the deconvolution process is 
to use a computer in place of a CSR spectrometer. By simulating a spectrum 
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Fig. 36. Energy in the radiation pulse detected within central cone. Curves 1, 2, 
and 3 correspond to a local energy spread of 5, 10, and 15 keV at the entrance to 
the bunch compressor (see Fig. 33). The bunch charge is 3 nC 

rather than using an actual CSR spectrometer to record it, we have a complete 
control over such factors as resolution, shape and noise. Also, we know what 
the perfectly deconvolved profile function would look like. In the deconvolution 
of actual CSR spectral data, the presence of noise is usually the main limiting 
factor: the effects of noise on deconvolution are demonstrated in the following. 
For the purpose of examining the deconvolution process, we begin with noisy 
data, which of course, can be realized in a simulation process. When other 
aspects of deconvolution, such as errors in the prior information are examined, 
noiseless data are used. 

Our goal is to discover the limiting sensitivity of the measurement technique 
considered. Attention is therefore turned to the crucial question of how accu- 
rately the bunch profile parameters can be measured by this technique when 
the errors in the input data cannot be neglected. Surely, the preferable way 
to answer the question about errors in bunch form-factor measurement is the 
computer experiments way. The effects of various errors are included in a full 
simulation of deconvolution procedure. One important question to answer is 
about how well the modulus of bunch form factor must be known in order to 
perform a successful deconvolution. To obtain the sensitivity to error effects, 
we use a set of input error data for a simulated bunch profile spectrum (see 
Fig. 37). A 40% rms error (variable over the undulator scan) is included. The 
sensitivities of the final bunch profile function to the simulation input errors 
are summarized in Fig. 38. Only one random seed has been presented here 
as an example, but several seeds have been run with similar success. Fig. 38 
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Fig. 37. Electron bunch form factor corresponding to the case of 10 keV energy 
spread at the entrance of the bunch compressor. The solid curve presents an ideal 
measurement. Symbols present five sets of measurements with relative error ±40% 
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Fig. 38. Current distribution reconstructed from five sets of measurements with rel- 
ative error ±40% (see Fig. 37). The solid line corresponds to an ideal measurement. 
The tail time constant in all cases is 9 ps 

shows that rms of CSR measurement should be within ±50%. It is seen that 
the constrained deconvolution method is insensitive to the errors in the bunch 
form-factor. 
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Fig. 39. Current distribution reconstructed from an ideal set of experimental data 
for the form-factor (see Fig. 37), but with different values of time constants, 7 and 
11 ps, for the bunch tail (dotted lines). The solid line corresponds to the correct 
value of the time constant, 9 ps 

As shown above, the bunch profile function is given by (47). An important 
issue to deal with consists in determining the dependence of the deconvolved 
function on the tail constant. One way to investigate this matter is to use a 
set of tail constants which differ from the real one in the deconvolution of a 
simulated spectrum. Fig. 39 shows the results of such a test. At first glance 
all three traces of Fig. 39 are practically identical, indicating that the tail 
parameter is less important than it would have probably been expected. We 
can therefore conclude that the width of the tail parameter should be within 
±50% of the proper value. However, if the correct tail function is not readily 
determined (by a streak camera) there is some justification for deconvolution 
with a "good guess" because even the extreme case shown in Fig. 39 is superior 
to trace the profile function without prior information. 



6 Concluding remarks 

Electron pulse-shape measurement is achieved by a three-step process. 

(i) In the first step spectrum of coherent radiation P(u) is measured. Coherent 
radiation from dipole magnet is emitted over a wide spectrum and should 
be analyzed with a spectrometer or interferometer. The proposed technique 
of using CSR from an undulator may be a more promising approach. The 
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basic measurement in this case is extremely simple: the energy of the coherent 
radiation pulse within the central cone is recorded by a bolometer. In order to 
use the radiation from the undulator it has to be filtered: for such a purpose, 
a simple pinhole can be used. 

(ii) The second step is to deduce the modulus of the bunch form-factor from 
the equation P(ou) = p(uj)N 2 | F(u) | 2 ; p(ou) should be applied to the result 
of the measurement, to obtain the bunch form factor. A careful analysis of 
the diagnostic instrumentation is required to find its exact p(ou) function. The 
computation of the single-particle spectrum for practical accelerator applica- 
tions is a rather difficult problem. The long wavelength synchrotron radiation 
from bending magnets in the near zone integrates over many different vacuum 
chamber pieces and the quantity p(u) is usually very difficult to know with 
good accuracy within the long wavelength range. In principle, modern com- 
puters allow one to perform direct calculation of the factor p(ou). The results 
of such simulations will depend on a large number of parameters. They will 
provide the possibility of obtaining a numerical answer for a specific set of 
input data, but hardly help to understand the physics of the problem. Success 
come to those who start from the physical point of view and begin by making 
the right kind of approximations, knowing what is large and what is small in 
a given complicated situation. We hope that the simple model employed here 
offers a useful introduction to the more formidable problems which arise in 
the computation of spontaneous emission spectrum distortions for real accel- 
erators. 

(iii) The third step consists in the reconstruction of the electron bunch profile 
from the measured data. The CSR spectrum measurement has an important 
limitation: it does not provide any information on the phase content of the 
pulse. It is clear that the retrievability of the information about the electron 
bunch profile function will in general not be complete, for it is the particular 
measured square modulus of the form-factor that can be obtained, not the 
complex form-factor itself. No spectral phase information is available, and the 
complete profile function cannot be obtained in general. Nevertheless, the less 
complete information can be exceedingly useful in many cases. For example, 
downstream of the bunch compressor the bunch charge distribution is strongly 
non-Gaussian with a narrow leading peak and a long tail. In this case detailed 
information about the longitudinal bunch distribution can be derived from the 
measured square modulus | F(u) | 2 by measuring the bunch tail constant with 
a streak camera and by using all a-priori available information about profile 
function. In this paper we demonstrate that constrained deconvolution is a 
useful experimental tool testable at each step in its development. Although 
the present work is concerned primarily with fully compressed electron bunch 
downstream of the first bunch compressor, its applicability is not restricted 
to that area only: for example, multistage bunch compression is a suitable 
candidate for treatment by the methods described here. 
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